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Abstract 

The operator product expansion of current correlators at short distances, and the notion of 
QCD-hadron duality are the cornerstone of QCD sum rules. The extension of this programme 
to T 7^ is discussed, together with applications to hot hadronic propagators. Indications are 
that the hadronic spectrum suffers a substantial rearrangement with increasing temperature, 
| and hint on the existence of a quark deconfming phase transition. Phenomenological order 

parameters to characterize this phase transition are discussed. 
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OS ' Some time ago Bochkarev and Shaposnikov Q proposed an extension of the QCD sum rule 

program to non-zero temperature, and made an application to the two-point function involving 
the vector current. This application was reconsidered in Later on we discussed the axial- 
vector channel [|| and the nucleon channel 0] using Finite Energy QCD sum rules (FESR). 
Results from these analyses indicate a substantial rearrangement of the hadronic spectrum 
with increasing temperature, and hint on the existence of a deconfming phase transition. This 
was later confirmed in []|], where a formalism valid even for T near the critical temperature 
^ ■ was used. An extension of QCD sum rules to finite temperature entails the assumptions that 

(a) the OPE continues to be valid, except that now the vacuum condensates will develop an 
(a-priori) unknown temperature dependence, and (b) the notion of QCD-hadron duality also 
remains valid. At T=0, the thermal behaviour of the vacuum condensates is not calculable 
analytically from first principles; some model or approximation must be invoked, e.g. the dilute 
pion gas approximation, lattice QCD, etc.. At T = 0, QCD-hadron duality is well established 
phenomenologically, while the validity of the OPE in QCD (in the presence of non-perturbative 
phenomena) is made plausible after studying other field theories which can be solved exactly 
||. At T ^ 0, it is important to examine critically the basis for the two assumptions above, 
especially in the light of some recent adverse criticisms 0. 

The basic object under consideration is the retarded (advanced) two-point function, involving a 
local current operator J(x), after appropriate Gibbs averaging 
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U(q,T)=i J (fix exp(iqx) 0(x o ) « [J(x), J f (0)] » , (1) 
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where 

<< A ■ B »= exp{-E n /T) (n\A ■ B\n) /Tr(exp{-H/T)) , (2) 

n 

and | n > is a complete set of eigenstates of the (QCD) Hamiltonian. The OPE of IT(q, T) is 
formally written as 

II (q, T) = d«I»+J2 C r (q) « O r » , (3) 

r 

The unit operator / in Eq.(3) represents the purely perturbative piece, and the OPE is assumed 
valid, even in the presence of non-perturbative effects, for q 2 < (spacelike), and | q 2 | 3> A.q CD . 
The states | n > entering Eq.(2) can be any complete set of states, e.g. hadronic states, quark- 
gluon basis, etc.. If one wishes to extend smoothly the T = QCD sum rule program to finite 
temperature, then the natural choice for this set is the quark-gluon basis, as first proposed in 
0]. I discuss below some arguments in support of the validity of the OPE at T / 0, Eq.(3). 
Assuming this, and invoking duality at T = 0, one now raises the temperature by an arbitrary 
small value. The hadronic spectrum and the expectation values of the QCD operators in the 
OPE will hardly change. Hence, it is reasonable to expect that the inter-relationship between 
QCD and hadronic parameters effected by duality will remain valid. An abrupt dissapearance 
of this inter-relationship as soon as T ^ 0, as advocated in JjJ, appears as an unlikely possibility. 
One should keep in mind that the current in Eq.(l) is an object external to the heat bath, and 
thus need not be in thermal equilibrium with the particles in the medium. The fact that at low T 
these particles are predominantly hadrons, is irrelevant for the argument. The external current 
will still convert to quark-antiquark pairs, which contribute e.g. to the perturbative operator / 
through loops, where they can have any value of momentum. In addition, there will be thermal 
non-perturbative effects as parametrized by the condensates. These are always distinguishable 
from perturbative effects, since they have a different q 2 dependence. In many cases they have 
also a different T-dependence, and may contribute only to the real part of the correlators. As 
the temperature is increased further, the hadronic spectrum will begin to change shape, in pace 
with the thermal behaviour of the expectation values of the various operators in the OPE. I shall 
discuss later how these changes emerge from QCD sum rules in agreement with expectations 
from general physical considerations. 

I turn now to the validity of the OPE at T / 0. Just as at T = 0, no rigorous proof can be 
given since one cannot solve QCD analytically and exactly. Instead, other field theory models 
which can be solved exactly have been used to argue for the validity of the OPE at T = |J. 
These models include the O(N) sigma model in the large ./V limit, and the Schwinger model, 
both in two dimensions. Comparing the short distance expansion of the exact solution for a 
Green function with the OPE counterpart, one finds that they are identical §. One can also 
show that this is also the case at finite temperature ||, as I outline next. Let us consider first 
the O(N) sigma model in 1+1 dimensions which is characterized by the Lagrangian 

C = \[d,a a {x)][d,a a {x)], (4) 

where a = 1,...N and a a a a = N/ f, with f being the coupling constant. In the large N limit this 
model can be solved exactly (for details see Q), it is known to be asymptotically free, and in 
spite of the absence of mass parameters in Eq.(4), it exhibits dynamical mass generation. In 
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addition, in this model there are vacuum condensates, e.g. to leading order in 1/N : (0|a|0) = 
y/N vn? , whereas all other condensates factorize, viz. (0|a fc |0) = (y/N m?) k . The a field is: 



a = f{d^a a ) 2 /\/N, and we are interested in the Green function associated with the propagation 
of quanta of this a field. We have calculated this Green function at finite temperature. Its 
imaginary part can be integrated analytically in closed form and is 

Imr(w,q = 0,T) = -L [l + 3n B (w/2T)] 



+ 



2uj 

1 
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« a » 6 << a 2 » 

+ T7 5 + 
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(5) 



where the first term above corresponds to the perturbative contribution, the second to the non- 
perturbative, and ns is the thermal Bose factor. Equation (5) is valid in the time-like region; 
the space-like region counterpart vanishes in 2 dimensions. Since the model is exactly solvable, 
the thermal behaviour of the vacuum condensates can also be calculated, viz. 

<< a »=< a > [l + 3n B (u/2T)\ (6) 

In this case the vacuum condensates contribute to the imaginary part, and as Eq.(5) shows, the 
thermal dependence of the perturbative piece cannot be absorbed into the condensates. Hence, 
no confusion should arise between perturbative and non-perturbative contributions. 
Finally, I consider the Schwinger model in 1+1 dimensions, with the Lagrangian 

C = -~F ltv + #i 7 „ if, (7) 

where T>a = idn + eA^. This model has been solved exactly, and in the framework of the OPE 
||. The short distance expansion of the exact solution coincides with that from the OPE |J. 
Here, we are interested in the two-point functions 

H ++ (x) = (0|T{j+(x)j + (0)}|0) (8) 

H + _(x) = (0|T{j+(x)r(0)}|0) (9) 

where the scalar currents are: j + = V'l , j~ = i>L V^R i with i^l,r = (1 ± 75)^/2. The 
function n ++ (Q) vanishes identically in perturbation theory, and the leading non-perturbative 
contribution involves a four-fermion vacuum condensate. We have calculated the thermal be- 
haviour of these current correlators and obtain, e.g. for their imaginary parts in the time-like 
region (again, there is no space-like contribution in 2 dimensions) 

Imn ++ (w,q = 0,T) = (10) 

Im n + _(w, q = 0, T) = - [1 - 2n F {uj/2T)} (11) 

Hence, the choice of the fermion basis in the Gibbs average of current correlators does not imply 
confusing these fermions with condensates, as argued in As Eqs.(10) and (11) indicate, 
(perturbative) fermion loop terms and (non-perturbative) vacuum condensates develop their 
own temperature dependence, which in this particular example happen to be different. 
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Turning to applications, an inspection of a typical hadronic spectral function at T = shows 
some resonance peaks at low energy, followed by a hadronic continuum starting at some threshold 
so- Because of asymptotic freedom, this continuum is well approximated by perturbative QCD. 
Raising the temperature, one would expect resonance melting. This may be accomplished by 
the imaginary part of hadronic propagators growing with temperature, as first proposed in |J. 
As the resonance peaks become broader, one would expect sq to decrease. Close to the critical 
temperature for deconfmement, the hadronic spectral function should be rather smooth, and 
represented almost entirely by the quark-gluon degrees of freedom over the whole energy range. 
An indication for some decrease of sq with increasing T was obtained from QCD sum rules for 
the vector channel M -H, |10|| . However, it is the axial channel the one that provides conclusive 
evidence for this behaviour [||, ||. In fact, the solution to the FESR gives 

-o(o) - mo) • 1 ) 

The temperature behaviour of the imaginary part of a hadronic propagator is difficult to obtain 
from QCD sum rules, since these are not too sensitive to resonance widths. Independent calcu- 



lations in other theoretical frameworks [11] show that the pion, the nucleon, and the rho-meson 
propagators have imaginary parts which grow monotonically with T. Over a certain range of 
temperatures these behave as 

r(T) = IT ^_ F . (13) 

Thus, the asymptotic freedom threshold so(T), and the resonance width T(T), are suitable 
phenomenological order parameters to characterize the quark deconfmement phase transition. 
It should be mentioned, in closing, that there are some unresolved problems in finite temperature 
QCD sum rules, e.g. when one tries to use them to extract the thermal behaviour of the vacuum 



condensates [1C]. An unwelcome implication is a breakdown of the FESR program beyond the 
lowest moment sum rule, and hence a breakdown of Laplace transform QCD sum rules. A 
satisfactory resolution of this problem is eagerly awaited. 
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